Search for the Skyrme-Hartree-Fock Solutions for Chiral Rotation in N=75 Isotones 
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A search for the self-consistent solutions for the chiral rotational bands in the N = 75 isotones, 
^■^"Cs, ^'^^La, ^'^''Pr and ^'^^Pm is performed within the Skyrme-Hartree-Fock cranking approach 
using SKM* and SLy4 parametrizations. The dependence of the solutions on the time-odd contri- 
butions in the energy functional is studied. From among the considered four isotones, self-consistent 
chiral solutions are obtained only in ^^^La. The microscopic calculations are compared with the ^^^La 
experimental data and with results of a classical model that contains all the mechanisms underlying 
the chirality of the collective rotational motion. Strong similarities between the HF and classical 
model results are found. The suggestion formulated earlier by the authors that the chiral rotation 
cannot exist below a certain critical frequency is further illustrated and discussed, together with 
the microscopic origin of a transition from the planar to chiral rotation in nuclei. We also formu- 
late the separability rule by which the Tilted- Axis-Cranking solutions can be inferred from three 
independent Principal- Axis-Cranking solutions corresponding to three different axes of rotation. 

PACS numbers: 21.30.Fe, 21.60.Ev, 21.60.Jz 



I. INTRODUCTION 

Since the original 1997 work of Frauendorf and Meng 
, the phenomenon of chiral rotation in atomic nuclei at- 
tracts quite a significant attention. The effect is expected 
to occur in nuclei having stable triaxial deformation, and 
in which there are a few high-j valence particles and a 
few high- J valence holes. The former drive the nucleus 
towards prolate, and the latter towards oblate shapes. 
The interplay of these opposite tendencies may favor a 
stable triaxial deformation. For such a shape, the valence 
particles and holes align their angular momenta along the 
short and long axes of the density distribution, respec- 
tively. However, the nuclear-bulk moment of inertia with 
respect to the medium axis is the largest, which favors 
collective rotation about that axis. Thus, the particle, 
hole, and collective angular momentum vectors are apla- 
nar, and may form either a left-handed or a right-handed 
set. In this way, the two enantiomeric forms may give rise 
to pairs of rotational bands, which are called chiral dou- 
blets. It is expected that the energy splitting between the 
partners in such doublets is very small, and in fact the 
authors of Ref. , who have analyzed the experimental 
results, have used essentially the argument 'by elimina- 
tion' - the bands were suggested to be chiral partners 
mainly because their properties could not be explained 
within other scenarios used by the authors. 

The first doublet band, later reinterpreted as chiral 1] , 
was found in 1996 by Petrache et al. in i34pj. gj jnjqw, 
about 15 candidate chiral doublet bands are known in the 
A « 130 region, and about 10 in the A « 100 region. The 
bands in the A Ki 130 nuclei are assigned to the simplest 
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chiral configuration, in which there is one proton parti- 
cle on the 7r/ixi/2 orbital, and one neutron hole on the 
1^^11/2 orbital. Configurations for A 100 nuclei usually 
involve one 7rgg/2 proton hole and one i^/iii/2 neutron 
particle orbitals. A few cases with more than one active 
particle or hole were also found 4]. So far, experimental 
information about absolute B(E2) and B(M1) values for 
transitions within the observed bands is available only 
for ^^*Cs and ^^^La, from recent lifetime measurements 
by Grodner 0,01 and Srebrny 7], and collaborators. 

On the theoretical side, chiral rotation has been exten- 
sively studied by using various versions of the Particle- 
Rotor Model (PRM), in which the nucleus is repre- 
sented by the valence particles and holes coupled via the 
quadrupole-quadrupole interaction to a rotator, often de- 
scribed within the Davydov-Fillipov model |^ with mo- 
ments of inertia given by the irrotational-flow formula , 
see, e.g., Refs. 0,0,^3- However, the main concept of 
rotational chirality summarized above, came from consid- 
erations within the Frauendorf 's mean-field Tilted- Axis- 
Cranking (TAG) model jlS], which is used in parallel 
with the PRM. That model is a straightforward general- 
ization of the standard cranking approach to situations 
where the axis of rotation does not coincide with any 
principal axis of the mass distribution. 

Up to noWjall the TAC calculations for chiral rotation 
(see Refs. |3, H 13 13 fo'' examples) were performed 
within a phenomenological approximation, in which the 
mean field is given by a simple model potential. A more 
fundamental description requires self-consistent methods, 
which would provide a strong test of the stability of the 
proposed chiral configurations with respect to the core 
degrees of freedom. Self-consistent methods are also nec- 
essary to take into account all kinds of polarization of 
the core by the valence particles and full minimization of 
the underlying energies with respect to all deformation 
degrees of freedom, including deformations of the cur- 
rent and spin distributions. Application of self-consistent 
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methods to the description of chiral rotation is the sub- 
ject of the present paper. 

Our study concerns four N = 75 isotones, ^^°Cs, ^^^La, 
^^"^Pr, and ^^^Pm, which are the first nuclei in which can- 
didate chiral bands were systematically studied ^J. We 
used the Hartree-Fock (HF) method with the Skyrme ef- 
fective interaction. The results were obtained for two 
Skyrme parameter sets, SLy4 [l3| and SKM* [l8|, and 
the role of terms in the mean field that are odd un- 
der the time reversal was examined. Calculations were 
carried out by u sing a new version of the code HFODD 
(v2.05c) [lEllflllllj which was specially constructed by 
the authors for the purpose of the present study. From 
among the considered four isotones, self-consistent chiral 
solutions were obtained in ^"^^La. A brief report on the 
results obtained in ^^^La was given in Ref. |22l|. 

The paper is organized as follows. In Section^we dis- 
cuss some characteristic aspects of the TAG calculations 
within the self-consistent framework. Section rill Al brieflv 
recalls previous studies on chiral rotation in the con- 
cerned isotones, and Section Fill Bl describes all technical 
details of the present calculations - in particular the way 
in which the role of time-odd nucleonic densities was ex- 
amined. Energy minima obtained for non-rotating states 
are listed in Section IIII CI In Section IIII Dl rotational 
properties of the valence nucleons and of the core are ex- 
amined within standard Principal- Axis Cranking (PAC). 
In Section llll El we solve a simple classical model of chiral 
rotation and show that such a rotation cannot exist below 
a certain critical angular frequency. The HF solutions for 
planar and chiral rotation are presented in Sections IIII Fl 
and IIII (II respectively. In Section Fill HI we demonstrate 
that our results obtained for the three-dimensional ro- 
tation can actually be represented as a sum of three in- 
dependent one-dimensional rotations about the principal 
axes. The values obtained for the critical frequency and 
the agreement of our results with experimental level en- 
ergies are discussed in Section FTVI In Appendix ^ we 
study response of the single-particle (s.p.) angular mo- 
menta to three-dimensional rotation, and introduce the 
notions of soft and stijf alignments. 



II. HARTREE-FOCK TILTED-AXIS-CRANKING 
CALCULATIONS 

So far, the TAG model has been described in the litera- 
ture only in its phenomenological variant Therefore, 
in this Section we give several details that are specific for 
its self-consistent implementation. The discussion con- 
cerns mainly the way of iteratively solving the HF equa- 
tions which is adopted in this work. 

As far as non-rotating states are concerned, the HF 
method consists in minimizing the expectation value of 
the many-body Hamiltonian, H = T+V, in the trial class 
of Slater determinants. Here, T is the kinetic-energy op- 
erator, and y is a two-body effective interaction. Equiv- 
alently, one can formulate the method in terms of the 



energy density functional E{p}, which is minimized with 
respect to the one-body density matrix p on which it 
depends, and this latter representation is used in the 
present study. Since H or the density functional E{p} are 
invariant under rotations in space, it is clear that the HF 
solution is defined only up to an arbitrary rotation. For 
each solution it is useful to introduce an intrinsic frame 
of reference, whose axes we define as principal axes of the 
tensor of the electric quadrupole moment of the mass dis- 
tribution. Due to the mentioned arbitrariness, this frame 
can be rotated with respect to the frame originally used 
to solve the HF equations, which we refer to as the pro- 
gram (or computer-code) frame. The program frame is 
the one defined by the axes x, y or z, used for solving the 
mean-field equations, e.g., in a computer code j3^. 

To describe rotational excitations within the TAG ap- 
proach, in the program frame one imposes a linear con- 
straint on angular momentum and minimizes the expec- 
tation value of the Routhian, 

H'^H-u}-l, (1) 

or the energy density in the rotating frame, 

E'{p} - E{p} - a; • Tr(/p) , (2) 

where / is the total angular momentum operator, and 
vector uj is composed of three Lagrange multipliers; it is 
called rotational frequency vector. Its components in the 
program frame are fixed as part of the definition of H' 
or E'{p}. Because of the rotational invariance, solutions 
obtained for the same length, but different directions of 
u: differ only by their orientation in the program frame, 
so that only the length, uj = \u:\, has physical meaning. 

Within the HF procedure, one obtains that the sought 
Slater determinant or the one-body density are built of 
the eigenstates of the s.p. Routhian, 

h' = h-u}-i , (3) 

where the mean-field Hamiltonian, h — T + T ^ 
SE{p}/Sp, is a sum of the kinetic energy, T, and the 
mean field, F. A phenomenological approximation to the 
HF approach consists in replacing F with a model poten- 
tial, like the Nilsson potential [23, whose deformation 
dependence is usually parametrized in terms of the mul- 
tipole deformations of the nuclear surface. Then, the 
expectation value of the Routhian, h' , is minimized over 
the deformation parameters of the potential within the 
standard Strutinsky method . 

Obviously, only the relative orientation of the angu- 
lar momentum vector with respect to the nuclear body 
carries physical information. Common orientation of the 
angular momentum vector and nucleonic densities with 
respect to the program frame can thus be arbitrary. In 
the phenomenological method, where the orientation in 
space of the nuclear surface is under control via the mul- 
tipole deformations ax^, one can take advantage of this 
fact and fix a2i = Im 0122 = so that the intrinsic and 
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program frames coincide. Minimization of the expecta- 
tion value of the Routhian, /i', at a given magnitude of uj 
is then performed by varying the direction of w and all 
except of 0:21 and Im a22- 

In the HF method, however, the Euler angles defin- 
ing the orientation of the intrinsic axes in the program 
frame are not free variational parameters, but compli- 
cated functions of densities, which in turn change from 
one HF iteration to another. The only possible way to 
make the two frames coincide is by imposing dynam- 
ical constraints on the off-diagonal components of the 
quadrupole tensor and requiring that they vanish. Then, 
one can vary the cranking frequency vector like in the 
phenomenological approach. This is actually the only 
way to proceed if the energy dependence on the intrinsic 
orientation of u) is sought. However, constraints strong 
enough to confine the nucleus easily lead to divergencies, 
and adjusting their strengths properly may be a cumber- 
some task. If only the energy minimum is of interest, a 
more natural and incomparably faster way is to fix uj in 
the program frame and let the mean potential reorient 
and conform to it self-consistently in the course of the 
iterations. The intrinsic axes now do become tilted with 
respect to the program frame. 

The Kerman-Onishi theorem [2^ states that in each 
self-consistent solution the total angular momentum vec- 
tor, I — (/), is parallel to a;. In calculations, the angle 
between those vectors converges to zero very slowly in 
terms of the HF iterations, because the whole nucleus 
must turn in the program frame in order to align its / 
with the fixed u). Therefore, a much faster procedure is 
to explicitly reset uj in each iteration to make it paral- 
lel to the current J, while keeping its length, cj, constant 
|20| . This purely heuristic procedure does not correspond 
to a minimization of any given Routhian. However, once 
the self-consistent solution is found, it is the Routhian 
for the final angular frequency that takes its minimum 
value. 

Some quantities, like mean angular momenta and mul- 
tipole moments, are easiest to discuss only when ex- 
pressed in the intrinsic frame of the nucleus, but it is a 
natural way to calculate them first in the program frame. 
Since the two frames do not necessarily coincide, one has 
to find the axes of the intrinsic frame by diagonalization 
of the quadrupole deformation tensor and to transform 
the considered quantities into that frame by use of the 
Wigner matrices. (Such a procedure may fail when the 
solutions have vanishing quadrupole moments, but in the 
present paper such cases are of no interest and will not 
be discussed). 

The HF TAG solutions are arbitrarily tilted with re- 
spect to the program frame, and, moreover, their orien- 
tation is not known a priori. Therefore, when solving 
the problem numerically one should ensure such condi- 
tions that the same solution be represented equally well 
in all orientations. In particular, the energy must not 
depend on the orientation. If the s.p. wave- functions are 
expanded onto a basis, this means that the choice of the 



s.p. basis and of the basis cut-off must not privilege any 
axis of the reference frame. In the case of the Cartesian 
harmonic-oscillator (HO) basis used in the present study, 
this by definition amounts to taking the three oscillator 
frequencies equal and including only the entire HO shells. 
To obtain a reasonable description of deformed nuclei in 
such non-deformed bases, the only way is to use suffi- 
ciently many HO shells. 

III. RESULTS 

A. Experiment and previous calculations 

It the present study, the HF solutions corresponding 
to chiral rotation were sought in four N = 75 isotones, 
i30Cs, "^La, i34pr, and ^^epni. These are the first nuclei 
in which candidate chiral bands were observed and sys- 
tematically studied ^ . For the most recent experimental 
data refer to [13 for ^^OCs, to 0, El 113 for ^^^La, to 
H m for i34pi.^ and to ^ for i^spm. Absolute values 
of spins and reduced transition probabilities were mea- 
sured only in ^'^^La. The possible chiral bands known in 
this nucleus are shown in Fig. ^ with full symbols. The 
lowest one (circles) is the yrast in positive parity. Two 
closely-lying excited bands, Bl and B3 (squares and dia- 
monds) are known 0,0,01 and it is unclear which one of 
them should be viewed as the chiral partner of the yrast 
band. The yrast and Bl bands are of positive-parity with 
known spins, while for band B3 the spin assignments are 
tentative. 

For all the TV = 75 isotones in question, PRM calcu- 
lations were performed 0, 0, |2£, ^QJ , and generally a 
good agreement with experiment could be obtained for 
both the energies and the B(M1)/B(E2) ratios, with a 
proper adjustment of the model parameters. Lifetime 
measurements in ^'^^La 0, S revealed, however, that 
absolute values of the reduced transition probabilities sig- 
nificantly deviate from these earlier PRM predictions. 

The phenomenological TAG calculations were carried 
out for all the considered nuclei as well 

i 111113 The 

quadrupole deformations e of 0.16, 0.175, 0.175, 0.195 
and triaxiahties 7 of 39°, 32°, 27°, « 26° were found, 
respectively for ^^OCs, ^^^La., ^^^Pi, "Spm. Chiral so- 
lutions were obtained in a limited range of angular fre- 
quency. The lower limits corresponded to the critical 
frequency, o-icriti discussed in the present work. Only 
Refs. [iJI and [13 quote their values, which are Wcrit = 
0.3MeV/;i for ^^'^Pr and Wcrit < 0.2MeV/;i for ^^^Pui. 
Comparison with experimental energies was given only 
for 434pp [1^3, where the behavior of the bands was 
reproduced, on the average. 

B. Details of the calculations 

The present study is limited to the simplest chiral s.p. 
configuration, i^h\^^2 ^^11/2' "^^ich has been assigned in 
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FIG. 1: (color online). Energies obtained from the HF TAC calculations, open circles and plus symbols, and from the classical 
model presented in Sec. HI.D, dashed and solid lines, compared with the experimental data on the candidate chiral partner 
bands in ^'^'^La, full symbols. The HF results obtained for the SLy4 and SkM* forces are shown for the A*', G, and T time-odd 
fields included (see text). 



the literature to the candidate chiral bands observed in 
the N — 75 isotones. In order to check whether the 
existence of the HF chiral solutions is not a particular 
feature of one Skyrme parameter set, all calculations were 
repeated with two parametrizations, SLy4 |0 and SKM* 
|l8|. The parity was kept as a conserved symmetry, and 
the pairing correlations were not included. 

In forming the chiral geometry, orientations of several 
angular-momentum vectors play a crucial role. Their be- 
havior depends, among others, on the nucleonic densities 
which are odd under time reversal, like the current and 
spin densities. Therefore, taking those densities and cor- 
responding terms in the energy density functional into 
account seems important for the microscopic description 
of the chiral rotation. Self-consistent methods are best 
suited for such a task, and investigating the role of the 
time-odd densities was one of our priorities in the present 
study. 

The Skyrme energy density functional depends on 
time-even and time-odd nucleonic densities with coupling 
constants Cf , Cf C; , C/, C^-^ (10 time-even terms) 



and with couphngconstants C|, , C{, C^^ (10 

time-odd terms) |31|- The index t — 0,1 denotes the 
isoscalar and isovector parts. In standard parametriza- 
tions, which are used in the present work, and C| 
additionally depend on the isoscalar particle density. If 
the assumption of the local gauge invariance is made |3lj| , 
there are several following relations between the time- 
even and time-odd coupling constants, 

ci=-c;, = c7' = +c7'. (4) 

In the present calculations, the coupling constants C/ 
of the time-even terms for t = 0,1 were always set to zero, 
like in the original fits of the forces SLy4 and SkM*. In 
order to conform to the local gauge invariance, whose 
consequence is Eq. (0J , the coupling constants of the 
time-odd terms were set to zero, too. Apart from C/, all 
other time-even coupling constants were taken as they 
come from the parameters of the Skyrme force. 

To examine the role of the time-odd densities, we have 
performed three variants of calculations, throughout the 
text denoted as N, G, and T, and defined in the following 
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way: 

• Variant TV: All time-odd coupling constants are set 
to zero and the local gauge invariance conditions in 
Eq. Q are violated. 

• Variant G: Coupling constants C/ and C^-'aie 
taken as they come from the parameters of the 
Skyrme force and as required by the local gauge in- 
variance all other time-odd coupling constants 
are set to zero. 

• Variant T: Apart from Cf, all other time-odd cou- 
pling constants are taken as they come from the 
parameters of the Skyrme force and as required by 
the local gauge invariance 

The density-dependent and independent components of 
were suppressed or not simultaneously. Setting or 
not some time-odd coupling constants to zero implies ex- 
cluding or including in the calculations the corresponding 
time-odd terms of the mean field. In variant N, the mean 
field contained only the time-even contributions (apart 
from the cranking term) as in the case of the phenomeno- 
logical mean potentials. 

The calculations were carried out by using the code 
HFODD (v2.05c) [HSIl^. The code expands the s.p. 
wave-functions onto the HO basis, and uses the iterative 
method to solve the HF equations. Twelve entire spher- 
ical HO shells were included in the basis. It has been 
verified that increasing this number up to 16, changes 
the quantities important for the present considerations 
(deformation, moments of inertia, alignments etc.) by 
less than 1%. 



C. Energy minima in the N = 75 isotones 

As the first step, the HF calculations without crank- 
ing were performed to find the Tr/ijj^yj ^^11/2 bandheads. 
Obviously, the energetically most favored state of this 
configuration is obtained if the valence proton particle 
and neutron hole occupy the lowest and highest levels of 
the hii/2 multiplets, respectively. Table Ogives the ob- 
tained (3 and 7 deformations for each isotone and Skyrme 
parameter set. They practically do not depend on the in- 
cluded time-odd terms, and remain almost constant when 
cranking is applied later on in our calculations. Note that 
the present values of (3 are up to 1.5 times larger than 
those used in the previous phcnomcnological TAG cal- 
culations quoted in Section IIII Al Also the values of 7 
are more distant from the maximum triaxiality of 30° as 
compared to the earlier results by other authors. 

In i34pj. g^jjjj ^^^Pm, two minima with the same 
^^11/2 ^^11/2 configuration were found, which differ by 
the occupation of positive-parity states. The energeti- 
cally lower minima have similar positive-parity s.p. struc- 
ture as in ^"^"Cs and ^"^^La, but they correspond to almost 
oblate shapes of 7 = 53°-to-58°. The other minima have 



7 = 19°-to-22°, thus corresponding to triaxial shapes. In 
the following, those two kinds of minima in i34pj. ^^-^^^ 
i36pj-|^ are conventionally referred to as oblate and tri- 
axial. Such a structure of minima and configurations 
appears for both interactions studied here, SkM* and 
SLy4, and the corresponding sets of results are very sim- 
ilar to one another. Therefore, to save space, below only 
the SkM* results are shown for the triaxial minima, and 
only the SLy4 results for the oblate ones. 

In Section ^ we introduced the intrinsic frame of a 
nucleus as formed by the principal axes of the quadrupole 
tensor. Below, by the short (s), medium (m), and long (/) 
axes of our triaxial solutions we understand the intrinsic 
axes indexed so that {xD < (x^J < (x^), where Xi is the 
Cartesian coordinate for axis i = s, or 

D. Properties of the /iii/2 valence nucleons and of 
the core 

We began our study of rotational properties by per- 
forming the standard PAC calculations, in which we 
examined rotations of the found, triaxial HF solutions 
about their short, medium, and long axes. Let us recall 
that for the PAC about the principal axis i = s,m,l, 
not only the cranking- frequency vector, a;, but all the re- 
sulting s.p. and total mean angular-momentum vectors, 
j and /, have non-zero components, Wi, ji, only on 
that axis. 

Figure 13 gives the s.p. Routhians obtained from the 
PAC calculations in ^'^^La with the SLy4 force, without 
the time-odd fields. Routhians for other isotones, forces, 
and time-odd terms included do not differ in their main 
features from the shown ones. The upper and lower parts 
show the proton and neutron levels, respectively, and 
the left, middle, and right panels correspond to crank- 
ing about the short, medium, and long axes. Positive- 
parity and negative-parity levels are marked with solid 
and dashed lines. The negative-parity levels belonging 
to the /111/2 multiplet are marked with crosses. The 
Routhians occupied by the valence /111/2 proton particle 
and neutron hole are marked with full and open circles, 
respectively. 

In both kinds of nucleons, the lowest levels of the ft.11/2 
multiplet split strongly for cranking about the short axis, 
and have almost zero slope for the two remaining axes. 
The highest /111/2 levels behave similarly, but split in the 
case of the rotation about the long axis. The intermediate 
levels split for cranking about each axis, but only weakly. 

For the PAC about the i-th axis, slopes of the s.p. 
Routhians, e', approximately translate into the s.p. 
angular-momentum alignments, ji, on that axis accord- 
ing to the well-known formula 



which holds almost exactly in the present case, because 
changes of the mean field with rotational frequency are 
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FIG. 2: Single-particle Routhians from the HF PAC calculations in La with the SLy4 force and no time-odd fields. Proton 
(upper half) and neutron (lower half) Routhians for cranking about the short, medium, and long axes are shown. Positive-parity 

and negative-parity levels arc marked with solid and dashed linos, respectively. Negative-parity levels belonging to the hii/2 
multiplet are marked with crosses. The Routhians occupied by the valence /111/2 proton particle and neutron hole are marked 
with full and open circles, correspondingly. 
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TABLE I: Quadmpole (3 and 7 deformation parameters, parameters of the classical model, Js,m,i [^^/MeV], Sg^i \h], classical 
estimates for the critical frequencies and spins, uif^ [MeV/7i] and I^'ff a-nd the full HF TAG results for those quantities, 
LO^Ji [MeV/R] and 7^?5| [Ti], for the Af = 75 i sotone s. The HF results with the SLy4 and SkM* forces are shown for the A^, G, 
and T variants of calculation defined in Sec. 1111 51 



nucleus 


force 


(3 7 




Ss Si 


clas rclas tit Tttt 
"-"crit -'crit "-"crit -'crit 


i^^Cs 


SLy4 N 
G 
T 


0.24 49 
0.24 49 
0.24 49 


4.81 29.2 17.0 
5.50 37.1 21.3 
4.16 29.4 19.9 


5.41 4.86 
5.45 5.16 
5.49 5.20 


0.46 12.8 
0.37 13.2 
0.59 16.8 


SkM* N 
G 
T 


0.23 47 
0.23 47 
0.23 47 


5.86 31.3 17.9 
6.55 36.7 21.1 
5.69 33.4 20.0 


5.43 5.01 
5.47 4.97 
5.49 5.14 


0.43 13.0 
0.37 13.0 
0.43 13.9 


^^^La 


SLy4 N 
G 
T 


0.26 46 
0.26 46 
0.26 46 


7.18 28.7 19.1 
8.45 36.0 23.7 
7.12 31.7 22.2 


5.44 4.90 
5.60 5.21 
5.64 5.26 


0.57 15.9 0.68 18.8 
0.47 16.4 0.60 20.3 
0.60 18.5 


SkM* N 
G 
T 


0.25 45 
0.25 45 
0.25 45 


8.19 30.8 20.3 
8.81 35.9 23.5 
8.37 34.0 22.4 


5.47 5.03 
5.60 5.06 
5.63 5.21 


0.54 16.0 0.62 17.8 
0.46 15.9 0.54 17.8 
0.50 16.5 0.58 18.5 


i34pj. 


SLy4 N 
oblate G 
T 


0.26 58 
0.26 58 
0.26 56 


6.11 25.4 
7.21 31.3 
3.68 29.8 


5.00 4.36 
4.92 4.84 
5.26 4.60 




SkM* N 
triaxial G 
T 


0.23 22 
0.23 22 
0.23 22 


18.5 28.1 20.7 
21.4 32.6 24.3 
20.7 30.8 24.8 


5.38 5.28 
5.46 5.44 
5.52 5.57 


0.91 25.0 
0.82 26.1 
1.08 32.7 


^^^Pm 


SLy4 N 
oblate G 
T 


0.25 53 
0.25 53 
0.25 52 


24.7 
30.3 
29.0 


4.62 
5.38 
5.05 




SkM* N 
triaxial G 
T 


0.22 19 
0.22 19 
0.22 19 


15.0 27.5 12.4 
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FIG. 3: Single-particle angular momentum alignments of the 
valence /iii/2 proton particle and neutron hole in ^"^^La, ob- 
tained from the PAG about the short, medium and long axes. 
The HF results with the SLy4 force are shown for the G vari- 
ant of calculation defined in Sec. IIII 51 

negligible. Thus, one can infer from Fig. [21 that the va- 
lence /111/2 proton particle and neutron hole have non- 
zero ahgnments, jf and j^", only for the PAC about short 
and long axes, respectively. This can be seen directly in 
Fig. 121 which shows jf and calculated as mean val- 
ues of the angular-momentum operator for the variant G 
of the calculations with the SLy4 force. Results for the 



short, medium, and long axes are plotted with dotted, 
solid, and dashed lines. 

It turns out that further important properties of the 
valence /111/2 particle and hole can be deduced from 
pure symmetry considerations, which we present in Ap- 
pendix ^ We consider an isolated two-fold degenerate 
s.p. level in a fixed potential symmetric with respect to 
the D2 group (rotations by 180° about the three princi- 
pal axes). The main conclusion relevant for the present 
case is the following. If in the PAC a s.p. state has a non- 
zero alignment, ji, only for rotation about one principal 
axis, then in the TAG its angular- momentum vector, j, 
will point along that axis and remain approximately con- 
stant regardless of the length and direction of the applied 
cranking frequency, u). We call it a stiff alignment. Thus, 
the angular momenta, and j'*, of the valence /iii/2 pro- 
ton particle and neutron hole are stiffly aligned along the 
short and long intrinsic axes, respectively, as expected 
for the chiral geometry. 

In our self-consistent calculations, the s.p. spectrum of 
the mean- field Hamiltonian, h of Eq. (jSJ, exhibits the 
two-fold Kramers degeneracy only in variant N of the 
calculations, when no time-odd fields are taken. It seems 
that the lowest and highest /iii/2 levels can be indeed 
treated as isolated: their angular-momentum coupling to 
other s.p. states is weak, which can be seen from the 
small curvature of their PAC Routhians in Fig. [3 The 



8 



assumption about fixed potential is also justified, because 
in our results the deformation remains nearly constant 
with rotational frequency. 

In the present case, the remnant coupling to other 
states and changes of the mean field alter the ideal pic- 
ture in that and ji'' are not strictly constant, but show 
some remnant dependence on the cranking frequency. As 
illustrated in Fig. 13 this dependence is to an approxima- 
tion linear, and the PAC alignments can be written in 
the form 

f ~ sj, + SJPu; , ~ siii + Sj'^u, , (6) 

where is and ii denote the unit vectors along the short 
and long axes, respectively. The quantities Sg and si are 
initial alignments at vanishing frequency, and 5J^^ and 
5J^ are tensor coefficients, representing the s.p. contri- 
butions to the total inertia tensor, , of the nucleus. 

The PAC allows to estimate the diagonal components 
of (5i/, which are slopes of the ji{uji) curves plotted in 
Fig. 13 One can see from the Figure that they may attain 
up to ~4?i^/MeV. This is a significant value compared 
to the total moments of inertia, discussed below and col- 
lected in Table For ^^'^La, the total moments are of 
the order of ~8-36 ?i^/MeV, depending on the axis. 

The same PAC calculations provide the total align- 
ments. Is, Im, h, on the short, medium, and long axes, 
respectively. They are plotted in Fig.^jwith dotted, solid, 
and dashed lines, correspondingly, for all the considered 
cases. In the nearly oblate minima in i34pj. g^^^i^ ^^^Pm, 
an attempt to crank around the medium axis leads to 
such a self-consistent readjustment of the shape that so- 
lutions corresponding to the rotation about the long axis 
are obtained. 

The bends in the curves in Fig. 01 like the ones for 
Is in ^'^°Cs, are caused by smooth crossings of the s.p. 
levels. Otherwise, the Ii{uji) dependence is linear to a 
good approximation, and therefore the rotation can be 
called rigid. The corresponding slopes give the collective 
total moments of inertia, J7s, J^m, and with respect 
to the short, medium, and long axes. They contain the 
valence particle and the valence hole contributions, 6J^p 
and 6j'^, defined in Eq. 

At zero frequency, the cranking around the medium 
axis gives a vanishing angular momentum, while the 
cranking around the other two axes give non-zero lim- 
iting values, equal to the initial alignments, Ss and s;, of 
the odd particle and hole. In view of the considerations 
presented later on in this paper, we have extracted the 
actual values of the parameters J's, Jm-, Ji-, and Ss, s;, 
by fitting straight lines to the calculated alignment vs. 
frequency curves shown in Fig. 0] Whenever there was 
a bend in the calculated dependence, the line was fitted 
in the rotational frequency range below the bend. Since 
the alignment on the short axis for the oblate minimum 
in ^36pni shows a particularly complicated behavior, we 
have not assigned any value to J7s in this case. The ob- 
tained values of the parameters are listed in Tabled and 
discussed in Section Hvl These results confirm that the 



moment of inertia with respect to the medium axis is the 
largest. 

The PAC method allows for estimating the diagonal 
components, J7s, J7m, and of the inertia tensor, J^. In 
order to examine, as far as possible, the off-diagonal com- 
ponents we performed a kind of perturbative test within 
the TAC method. We applied to the non-rotating self- 
consistent solutions the cranking frequency vector, u), in 
several directions, performed only one diagonalization of 
the resulting s.p. Routhian ||2Jl, and investigated the re- 
sponse of the mean angular momentum, I. It turned 
out that the off-diagonal components of J' are negligibly 
small in all the considered isotones. 

The microscopic results presented so far suggest that 
the considered system can be modeled by two gyroscopes 
of spins Ss and s; rigidly fixed along the short and long 
axes of a triaxial rigid rotor characterized by the mo- 
ments of inertia, J7s, J^m, Ji, of which is the largest. 
It is instructive to solve the associated problem of mo- 
tion in the classical framework, which is done in the next 
Section. 



E. Classical model 

The classical model of chiral rotation was briefly in- 
troduced in Ref. [l^l; here we give its more detailed de- 
scription and discussion. In order to define the model, 
we begin with elementary considerations related to the 
dynamics of rigid bodies. By a classical gyroscope we 
understand an axial-shape rigid body with moments of 
inertia, J\\ and J71, with respect to the symmetry axis 
and an axis perpendicular to it, respectively. Such a body 
spins with fixed angular frequency f2 around its symme- 
try axis; one can imagine that this motion is ensured by 
a motor and frequency regulator such that 51 is strictly 
constant in time. 

Furthermore, let us imagine that the spinning body is 
rigidly mounted on another rigid body that has triaxial 
inertia tensor and three moments of inertia Jl, J!^, and 
with respect to its short, medium, and long axis, re- 
spectively. Then, the angular frequency fl is maintained 
fixed with respect to the triaxial body, irrespective of how 
the whole device moves, and the angular frequency vector 
O has by definition three time-independent components 
57s, 51m, and 51/. To simplify our considerations let us 
assume that the axis of the gyroscope coincides with the 
short axis of the triaxial body, i.e., 51<,=51, 51^=0, and 
51; =0. In this case the principal axes of the tensor of 
inertia of the whole device coincide with those of the tri- 
axial body, and the three moments of inertia of the device 
read, 

Js = J's + J\\ 1 
Jm — J'm^ J^, (7) 

Jl = Jl + J^. 

We assume that the device rotates with the total an- 
gular frequency vector a;, which in the moving frame of 
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FIG. 4: Total angular momentum alignments obtained from the PAC about the short, medium and long axes in "°Cs, "^La, 
"*Pr, and "^ Pm. T he HF results with the SLy4 and SkM* forces are shown for the A'', G, and T variants of calculation 
defined in Sec HTTSl 



the triaxial body has components w^, w™, and w;. In 
general, these components may vary with time, although 
later we study only such a motion of the device when they 
are time- independent. The kinetic energy of the device 



is the sum of that of the triaxial body and gyroscope. 



— ^triax ~t- ^gyro5 (8) 
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where 

Ttriax — 5 {Jl^"^ + J'n^m + ■^I'^'l) J (9) 

Tgyro - ^{j\\iiJs + ny+J^u;l + J^u;f), (10) 

hence 

We see that the total kinetic energy is a sum of three 
terms. The first one represents the rotation of the entire 
device irrespective of the fact that it contains a spinning 
gyroscope; it depends only on the total moments of in- 
ertia. The second one represents the additional energy 
coming from the spinning gyroscope and depends on its 
time-independent spin Sg = J7\\i^, while the third one is 
a constant which can be dropped from further consider- 
ations. 

It is obvious that if we add two other gyroscopes 
aligned with the medium and long axes and spinning 
with spins s,„ and si, respectively, the second term can 
be simply written as a scalar product ■ s, where vector 
s has components Ss, Sm, and si. In this case, the total 
moment of inertia will be a sum of contributions from 
the triaxial body and three gyroscopes. We note in pass- 
ing that exactly the same result is obtained for a single 
spherical gyroscope, jTji = J7j_, tilted with respect to the 
principal axes of the triaxial body in such a way that it 
has spin components equal to s^, Sm, and s;. 

The total angular momentum, J, of the system reads 

I = Juj + s, (12) 

where s is the above vector sum of spins of all the gy- 
roscopes and J^u) stands for the tensor product of the 
moment of inertia tensor with the angular frequency 
vector w. In absence of potential interactions, the La- 
grangian of the system is equal to the total kinetic en- 
ergy (|TT|l . generahzed to the case of three gyroscopes, 
and thus is given by the formula 

L ^ ^ivjuj + u} ■ s , (13) 

where we dropped all constant terms. Taking the lab- 
oratory components of a; as generalized velocities, it is 
easy to check that the generalized momenta are equal to 
the laboratory components of I. This fact allows us to 
write the Legendre transformation js^l and to obtain the 
Hamiltonian of the system, 

H = ujI -L = ^ujjuj . (14) 

Since the Lagrangian l|13|l docs not depend explicitly on 
time, the Hamilton function H14I) is a constant of mo- 
tion, and is identified with the total energy, E, of the 
system. Consider now a particular type of the Routhian 
[33, H', namely such that no variables undergo the Leg- 
endre transformation. In such a case, H' = —L, and by 
rewriting the Routhian in terms of the Hamiltonian one 
obtains 

H'^H-u}-I. (15) 



Equations of motion for the model can be derived in 
the following way. As for any vector, the time deriva- 
tives of the angular momentum vector, dtl - taken in 
the laboratory frame and d'^I - taken in a frame rotat- 
ing with angular frequency u>, are related by the formula 
dtl = d^I + u! X I .32j. Since the angular momentum 
is conserved in the laboratory frame, dtl — 0, one ob- 
tains the Euler equations 32] for the time evolution of 
the angular-momentum vector in the body-fixed frame, 

d'^I =-u}xI . (16) 

The mean-field cranking approximation can only ac- 
count for the so-called uniform rotations, in which the 
mean angular-momentum vector is constant in the in- 
trinsic frame of the nucleus, 9^ J = 0. Because of that, 
we restrict the classical model studied here to such uni- 
form rotations. Due to Eq. (|12|l . for uniform rotations 
also the angular frequency vector is constant in the in- 
trinsic frame. The Euler equations H16() now take the 
form a; X / = 0, and require that uj and / be paral- 
lel. The same condition holds for the HE solutions and is 
known as the Kerman-Onishi theorem; see SectionllTland 
[2^ . The Euler equations can now be easily solved for 
the considered classical model. However, to show further 
analogies with the HE method, in what follows we find 
the uniform solutions by employing a variational princi- 
ple. 

According to the Hamilton's principle, motion of a me- 
chanical system can be found by making the action in- 
tegral, J L dt, stationary. The real uniform rotations 
obviously belong to a wider class of trial motions with 
I and u} being constant in the intrinsic frame, but not 
necessarily parallel to one another. Within this class, La- 
grangian H13(l does not depend on time. Therefore, ex- 
tremizing the action for a given value of u;=|ci;| reduces 
to finding extrema of the Lagrangian as a function of 
the intrinsic-frame components of uj. Since H' = —L, 
the Routhian l|15|l can be equally well used for this pur- 
pose. This provides us with a bridge between the classical 
model and the quantum cranking theory, where an anal- 
ogous Routhian Q is minimized in the space of the trial 
wave- functions. 

Extrema of H' with respect to the intrinsic-frame com- 
ponents of a; at a given length of lu can be found by us- 
ing a Lagrange multiplier, — ^/i, for cj^ (the factor — i 
being added for later convenience). We continue further 
derivations for the case of two gyroscopes aligned along 
the s and I dictated by the microscopic results 

presented in Section Fill Dl i.e., we employ the classical 
model for Sm=0. Setting to zero the derivatives of the 
quantity 

+ i^i - Jl)ujf] - {uJsSs + t^isi) (17) 
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with respect to Ug, ujm, t^i, one obtains 

= Ss/{il-Js) , (18) 
W™(m- Jm) = , (19) 

= sl/{^Ji-Jl) . (20) 

Equation (|19l) gives either uj„i — or ji ~ Jm , leading to 
two distinct classes of solutions. 

If LOrn — then both a; and I lie in the s-l plane. This 
gives planar solutions, for which the chiral symmetry is 
not broken. All values of /i are allowed, and the Lagrange 
multiplier must be determined from the length of a; cal- 
culated in the obvious way from l(TS)l and (^0)1 . Figure Et- 
shows uj versus for sample model parameters, extracted 
from the ^^^La HF PAC solutions with the SLy4 force 
with no time-odd fields, and listed in Table HI The solu- 
tions marked as A and D exist for all values of w, while 
above some threshold frequency, wthr , two more solutions 
appear, B and C. This threshold frequency can be deter- 
mined by finding the minimum of io in function of /i, and 
reads 



Wthr 



\Jl-Js\ 



3/2 



(21) 



The value of wthr coming from the present HF calcula- 
tions is rather high, higher than 1 MeV/?i. Since bands B 
and C are situated far above the yrast line (see Fig. 01) 
they will not be subject of further analysis. 

For [i=Jm^ all values of are allowed, while com- 
ponents in the s-l plane are fixed at lOs — Ss/{Jm — Js) 
and oJi = Si/{Jrn — Ji)- Consequently, the angular mo- 
mentum has non-zero components along all three axes, 
and the chiral symmetry is broken. For each value of 
w, there are two cases differing by the sign of w^, and 
thus giving the chiral doublet. The fact that LOg and 
are constant leads to the principal conclusion that chi- 
ral solutions cannot exist for w smaller than the critical 
frequency 



clas 



\Jrri \J s 



Sl 



Jm — Jl 



-I 1/2 



(22) 



At that frequency, and with uj„i~Q, the chiral solution 
coincides with the planar band A. 

In the uj space, the four planar solutions form a hy- 
perbola in the s-l plane, while the chiral doublet cor- 
responds to a straight line perpendicular to that plane. 
These curves are shown in Fig. ^p. Figure gives the 
angular momentum in function of rotational frequency 
for all the presented bands. With increasing lo, the so- 
called dynamical moment, = dl/dw, asymptotically 
approaches Ji for bands A and B, and Jg for bands C 
and D. For the chiral band, / is exactly proportional to 
UJ with the coefficient J^,n. Thus, the critical spin, I[ 
corresponding to the critical frequency H22|l reads 



;las 



rclas 



J clas 
m'^crit 



Figure [Sji summarizes the energies in function of spin. 
The spin quantum number, J, is related to the length, 
/, of the angular momentum vector, J, by the condition 
/(I -I- 1) = |Jp. At low angular momenta, the yrast 
line coincides with the planar band D. Then it continues 
along the planar solution A. Since the moment of inertia 
J^m is the largest, beyond the critical frequency the chiral 
solution becomes yrast, thereby yielding good prospects 
for experimental observation. 

Altogether, the classical model described here is de- 
fined by five parameters, Js-, Jm, Ji, Sg, and si, which 
are extracted form the microscopic HF PAC calculations 
and listed in Table The model can then be applied to 
predict properties of the planar and chiral TAC bands, 
and these predictions can be compared with the HF TAC 
results. Such a comparison is presented and discussed in 
the following Sections. 



F. Planar solutions 

We began our self-consistent TAC calculations by find- 
ing planar solutions corresponding to the classical band 
A. The first point of each band of this kind was obtained 
by restarting the HF iterations from the previously con- 
verged non-rotating solution and by applying the initial 
cranking frequency vector with non-zero components on 
its short and long intrinsic axes. Once convergence was 
achieved, the obtained solution served in turn as the 
starting point for the next value of the rotational fre- 
quency. We proceeded in this way with the frequency 
step of 0.05MeV/?i. We followed each band diabatically, 
i.e., by exciting particles near the Fermi level whenever 
an empty and an occupied s.p. level of the same parity 
were about to cross, so that always the states with the 
same physical properties were occupied. 

In the solutions corresponding to a planar rotation, 
the cranking frequency vector had non-zero components 
on the short and long intrinsic axes in the self-consistent 
results. To give some insight into the s.p. structure of 
those solutions, in Fig. |^ we give the proton and neutron 
s.p. Routhians for the planar band in ^^^La obtained with 
the SkM* force and no time-odd fields. Parts of the plots 
to the left of the thin vertical lines are relevant for the 
planar results. Contrary to the PAC Routhians of Fig. |21 
now both the lowest and the highest /iii/2 levels split with 
rotational frequency. This is consistent with the picture 
that the valence particle and hole angular momenta, 
and j'^ , aligned on the short and long axes, respectively, 
now both have non-zero projections on the tilted axis of 
rotation. 

In order to examine the angular momenta of the va- 
lence nucleons, and j'*, in Fig. [7|we plot their projec- 
tions onto the short (dotted line), medium (solid line) and 
long (dashed line) intrinsic axes, for all the self-consistent 
solutions in ^'^^La. The positive and negative alignments 
are those of the lowest proton and highest neutron /111/2 



(23) levels; the latter can be considered as representing ~j 
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FIG. 5: (color online). Planar bands A, B, C, D, and the chiral band obtained from the classical model, a) Rotational frequency, 
Lj(^). b) Intrinsic-frame trajectory of the rotational frequency vector, uj. The chiral solution corresponds to a straight line 
perpendicular to the figure plane, and intersecting it at the marked point. That perpendicular direction represents u>m- c) 
Angular momentum, I{lu). d) Energy, E{I). 



Parts of the plots to the left of the thin vertical lines 
concern the planar bands. It can be seen that, indeed, 
the proton particle and the neutron hole align their an- 
gular momenta on the short and long axes. Further- 
more, those alignments change rather weakly with ro- 
tational frequency which means that the wave-functions 
are strongly confined by deformation (deformation align- 
ment). 

It is worth emphasizing that the intrinsic-frame trajec- 
tories of us along the self-consistent bands almost exactly 
follow the classical ones, in all the considered cases. This 
is illustrated in Fig. |S| where the dashed lines represent 
classical A bands with parameters of Table HJ and the 
HF results are marked with open circles. For the oblate 
bands in ^^^Pm the classical line is not shown because 
the parameters JT^ and Ss could not be unambiguously 
extracted from the PAC calculations; see Section Fill Dl 
However, the HF results follow a curve that very much 
resembles the classical hyperbola. For the triaxial min- 
ima in i34pj. g^jjj i36pj^^ self-consistent planar bands 
could be obtained because of multiple level crossings. 

Energy in function of spin also shows a striking agree- 
ment between the classical and self-consistent results for 
the planar bands. This can be traced in the case of ^^^La 
in Fig.^ by following the same symbols as those in Fig.|Sl 



Some deviations are visible only for rather high angular 
momenta. 



G. Chiral solutions 

The planar HF solutions were easily obtained by ap- 
plying small cranking frequency increments to the non- 
rotating state. For chiral bands, analogous task was 
more difficult, because these bands start at finite fre- 
quencies, which in the present case are not lower than 
w 0.5MeV/?i. Several level crossings may occur between 
uj = and such a high frequency, and it is difficult to 
spot the required s.p. configuration at high frequencies. 
A hint on how to follow the ^^11/2 configuration 

diabatically comes from the classical prediction that the 
chiral band branches off from the planar solution (at the 
point corresponding to the critical frequency). One can 
thus restart iterations from the planar band by apply- 
ing cranking frequency with non-zero component on the 
medium axis. 

As the first step we performed a kind of perturbative 
search along the planar bands, which turned out to be a 
very reliable test of where the chiral solution may exist. 
Such a search gives us also some understanding of why 



13 



Protons Neutrons 




Rotational frequency [MeV/^] 



FIG. 6: Similar as in Fig.|5|but for the proton and neutron s.p. Routhians from the HP TAG calculations with the SkM* force. 
The thin vertical line is drawn at LO^Jt- The Routhians to the left and to the right of this line correspond to the planar and 
chiral bands, respectively. 



the chiral solutions do not appear in several cases. The 
procedure we applied was the following. To each con- 
verged point of the planar band, a small additional com- 
ponent, Wm, of the angular frequency along the medium 
axis was added, with ujs and uji unchanged. The resulting 
s.p. Routhian ||2J) was diagonalized only once. Then, it 
was checked whether in the resulting non-selfconsistent 
state the angular momentum and rotational frequency 
vectors were parallel, as required by the Kerman-Onishi 
necessary condition of self-consistency; see Section^ We 
can guess that in nuclei stiff against deformation changes, 
the direction of I is the only degree of freedom, and thus 
the Kerman-Onishi condition is also sufhcient. If / is 
parallel to w in the non-selfconsistent state after one di- 
agonalization, then it is very probable that further iter- 
ations may lead to a converged chiral solution. Indeed, 
this was always the case, and never a chiral solution was 
obtained, in spite of several attempts, if that simple test 
gave negative result. 

The condition for / and uj being parallel can be written 
in the form 

I^^Irn^Il^ (24) 

Note that the Is/^s and Ii/uii ratios must be very close 
to each other in the non-selfconsistent state, because 



the Kerman-Onishi condition is fulfilled for the self- 
consistent planar state. Therefore, the test consists in 
checking for each point of the planar band if /„j/w,„ is 
equal to Is/i^s ~ Ii/uJi- In fact, this is reliable only if the 
time-odd fields are switched off, because in their pres- 
ence, the Im/^m ratio calculated perturbatively is sig- 
nificantly smaller than the self-consistent result would 
be. The reason is that the relevant components of the 
time-odd fields become active only after self-consistency 
is achieved for a non-zero LOrn- The test was made with 
bJm — 0.05 MeV/?i. The discussed ratios, calculated for 
all the HF planar bands found in the N = 75 isotones, are 
plotted in Fig. |^ Plus symbols and open circles denote 
the HF values of I„i/ujm and Is/^s ~ h/^i, respectively. 
To guide the eye, in the same Figure we also plotted the 
ratios Im/'^m and Is/tOs = Ii/uJi corresponding to the 
classical-model chiral and planar bands, respectively. 

It can be a priori expected that chiral solutions do not 
appear in the oblate minima in i^^pj. g^j^j ^^^Pm, because 
of insufficient triaxiality. Indeed, the calculated values of 
the Im/^m and Ig/^s ~ Ii/^i ratios exhibit a compli- 
cated behavior, and do not become equal to one another 
at any point. In "'^^'^Cs, as well as in the triaxial minimum 
in ^^^Pm, the two ratios clearly approach each other. It 
seems that the only reason why they do not attain equal- 
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FIG. 7: Angular-momentum alignments of the lowest proton (positive values) and highest neutron (negative values) /iii/2 levels 
on the short, medium and long intrinsic axes from the HF TAG calculations in ^'^'^La. The thin vertical line is drawn at iJ^J^x- 
The curves to the left and to the right of this line correspond to the planar and chiral bands, respectively. The HF results with 
the SLy4 and SkM* forces are shown for the A*', G, and T variants of calculation defined in Sec. If If 51 



ity is that the planar bands were not found up to suffi- 
ciently high frequencies, because of level crossings. Note, 
however, that the moment of inertia associated with the 
medium axis, Im/'^rm significantly drops with angular 
frequency, which takes Im/oJm away from Is/ujs ~ h/^i, 
and defers their equalization to higher frequencies. This 
effect is much weaker in ^'^^La, where the ratios do be- 
come equal, slightly above the point expected from the 
classical model. Indeed, self-consistent chiral solutions 
were found in this described below. 

After the first diagonalization of the perturbative test, 
the HF iterations were continued in each case to achieve 
self-consistency. In ^^°Cs, ^"^^Pr, and ^"^^Pm, as well as 
for low rotational frequencies in '^'^^La, the iterations con- 
verged back to planar solutions. The same result was ob- 
tained for different initial orientations of u> with respect 
to the intrinsic frame. This provides a strong argument 
that, for the concerned configuration, no self-consistent 
chiral solutions exist at low frequencies. In ^^^La, for 
Lj high enough, converged solutions were obtained with 
/ having non-zero components on all the three intrin- 
sic axes, which corresponds to chiral rotation. To ex- 
amine the chiral solutions independently of the planar 
ones, the found fragments of chiral bands were used as 
starting points to obtain solutions for lower and higher 



frequencies. Calculations were performed with oj step of 
0.02 MeV/?i. At a certain value of decreasing w, the pla- 
nar orientation of I was regained in the intrinsic frame, 
and the solution merged into the previously found pla- 
nar one. In a natural way, that junction value of w can 
be regarded as the Skyrme-HF result for the critical fre- 
quency, and is denoted in the following as lu^J^. Values 
of obtained in the present calculations are collected 
in Tabled and discussed in Section Hvl On the side of 
highest frequencies, chiral solutions were obtained up to 
a certain value of and all attempts to go higher caused 
the iterations to fall into a different minimum. This is 
probably due to multiple smooth crossings of occupied 
and empty levels, particularly in neutrons. 

The HF results presented so far corroborate the main 
prediction of the classical model, that chiral solution ex- 
ists only above a certain critical frequency, at which it 
branches off from the planar one. Also the intrinsic- 
frame trajectory of a; along the HF chiral band is almost 
a straight line parallel to the medium axis, as in the clas- 
sical case. This is demonstrated in Fig. ^1 which shows 
projections of u> on the s-m and l-m intrinsic planes for 
the HF (plus symbols) and classical (line) results. The 
only difference is that is a bit higher than , and 
the HF line is shifted along the planar band to higher fre- 
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quencies. Although the chiral solutions have been found 
in a rather narrow uj interval, of about 0.1MeV/?i, the 
accompanying increase in Wm is significant, from zero to 
about 0.4MeV/?i. This is so because ojs and coi are rela- 



tively large and almost constant along the chiral solution. 

Figure ini shows the proton and neutron s.p. Routhians 
for the chiral solution obtained in ^^^La, with the SkM* 
force and no time-odd fields. A thin vertical line is drawn 
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at the value of wJrit- The Routhians to the left and to 
the right of this line correspond to the planar and chi- 
ral bands, respectively. Note first that the planar and 
chiral Routhians do indeed coincide at oj^J^- The chiral 
Routhians do not seem to exhibit any particular behav- 
ior. At high frequencies, the Routhian occupied by the 
neutron hole enters into a region of high level density and 
the corresponding s.p. state mixes with other negative- 
parity states. Thus, it is doubtful whether the valence 
neutron hole can be identified with a single state, and 
we do not examine its s.p. properties. The marking in 
open circles is only tentative. However, the /iii/2 proton 
particle is still well separated. 



The alignments of the angular momentum, , of the 
^11/2 proton particle on the short, medium, and long 
intrinsic axes for the SLy4 and SkM* forces and N, G, 
T time-odd fields are shown in Fig. [7| As in Fig. El the 
vertical line separates the planar and chiral bands. The 
plot confirms the stiff character of those alignments in the 
chiral solutions. Since the ujg and uji components of the 
cranking frequency vector hardly change along the chiral 
band, also the considered alignments on those axes, 
and , are nearly constant. Only the projection on the 
medium axis, j^, increases due to the increase in LOm 
from zero to about 0.4MeV/?i. 



H. Separability of the TAG rotation 

To complete the presentation of our results, we here 
formulate and discuss the separability rule, by which our 
planar (2D) and chiral (3D) HF TAG solutions turn out 
to be simple superpositions of independent HF PAC (ID) 
rotations about 2 or 3 principal axes, respectively. 

Take, for instance, the total alignments on the princi- 
pal axes, Ii{ijJi)^ where i = s,m, I, as functions of the cor- 
responding components of the angular frequency. Such 
quantities were examined in Section 1111 Dl for the PAC 
calculations, but they can be equally well extracted for 
the planar and chiral solutions. A comparison of the 
PAC and TAC results in ^^^La. is shown in the lower 
part of Fig. ^] for the T variant of calculations with the 
SkM* force. The PAC values of Is, Im, and Ii are plot- 
ted with lines, respectively: dotted, solid, and dashed. 
The TAC results are marked with symbols, respectively: 
down-triangles, diamonds, and up-triangles. The left and 
right parts of the Figure show the TAC values of 
for the planar and chiral solutions, respectively, ft is 
obvious that the results obtained for the rotation about 
tilted axes are almost identical to those obtained from 
the independent PAC calculations. 

Unlike the angular momentum, the average value of the 
total TAC Routhian, {H') of Eq. I^), cannot be trivially 
decomposed into contributions from rotations about the 
three principal axes. Yet, for each TAC solution, char- 
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acterized by an angular frequency uj, we have its three 
intrinsic components, uji, and we can consider a sum, 
H', of the corresponding values of the PAC Routhians, 

H\lo)= J2 (AOM, (25) 

which we call equivalent Routhian. We then compare the 
difference relative to its value at uj=Q, 

AH'{lu) = H'{uj) - H'{0), (26) 

with the analogous difference, 

AH'{co)^{H')iu;)-{H')iO), (27) 

computed for the full TAG Routhian. In the upper part 
of Fig. ^ the differences and l(77|) are plotted in 
function of uj as lines and points, respectively. 

For the chiral solution, as {H'){0) we take the same 
value as for the planar case, because the chiral band can 
be regarded as a continuation of the planar one, as dis- 
cussed in Section llll Gl The equivalent Routhian can be 
plotted only in such a frequency range in which the PAC 
solutions are obtained for the corresponding components 



uji. In that range, AH' deviates from AH' by not more 
than 30keV, as it can be seen from the Figure, whereas 
the total Routhian itself drops by about 8 MeV between 
Lu = and 0.6MeV/?i. Therefore, we conclude that the 
equivalent Routhian, constructed out of the PAC solu- 
tions, reproduces the full TAC Routhian to a very high 
accuracy. 

Another piece of information that is contained in the 
TAC results, and not directly in the PAC results, is how 
the total angular frequency, w, should be distributed over 
the three principal axes for each given point of a band 
representing the rotation about a tilted axis. That, how- 
ever, can also be determined from the PAC solutions by 
means of the Kerman-Onishi requirement (see Sec- 
tionm that the angular-momentum vector must be par- 
allel to the angular-frequency vector. Indeed, this condi- 
tion is equivalent to the system of three non-linear equa- 
tions, 

IsiuJs)/(^s = M, 
/m(Wm)Mn = M, (28) 
Il{uJl)/uJi = ^1. 

By plotting the three PAC functions Ii{uji)/LOi on the 
same plot as functions of oji we obtain a planar or chi- 
ral solution whenever two or three, respectively, of these 
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functions cross a horizontal line. Then, values of Ui cor- 
responding to the crossing points define the requested 
distribution of the components within the total angu- 
lar frequency ui. The agreement of the PAC and TAG 
alignments shown in Fig. ^2 guarantees that the above 
procedure gives a correct result. 

We have thus demonstrated that within our HF re- 
sults, the general 2D or 3D rotation separates into two 
or three, in a sense independent, ID rotations about the 
principal axes. This is possible mainly because our solu- 
tions are very stiff against deformation changes with ro- 
tational frequency. Had they been soft, rotations about 
different (principal or tilted) axes could cause different 
shape polarizations that would prevent such a simple su- 
perposition of motions. Yet, whenever the motions are 
separable, the PAC calculations supplemented with the 
Kerman-Onishi condition are actually sufficient to de- 



scribe the 2D or 3D rotations. Although this rule has to 
be, in principle, confirmed numerically in each particu- 
lar case, it is plausible that it will hold in all analogous 
cases of stiff alignments, by which the difficult TAC cal- 
culations can be replaced by much more easily performed 
PAC calculations. 



IV. DISCUSSION 

The main point left for the discussion concerns the val- 
ues of the critical frequency: where do they come from, 
whether they depend on the Skyrme force used, how they 
are altered by the different time-odd terms of the func- 
tional, how they would be influenced by the inclusion of 
pairing, and what their relation to the experimentally 
observed bands is. In this Section we give some remarks 
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on these and related topics. 

Table U summarizes the values of the classical-model 
parameters, J7s, Jrm Ju Sg, and s;, extracted from the 
HF PAC results in Section Fill Dl of the classical critical 
frequency, Wcrft , calculated from Eq. (|^ . and of the crit- 
ical frequency, lo^^^, obtained from the HF TAG calcula- 
tions as defined in Section Fill Gl In addition, the Table 
gives the corresponding values of critical spins, if^f^ and 
^cr^' respectively. 

First of all, it can be seen from the examples of ^^°Cs 
and ^■^^La that the SLy4 and SkM* forces give quite sim- 
ilar values for all the concerned quantities. The differ- 
ences are not larger than variations within one force due 
to taking into account different time-odd fields. For all 
the four isotones in question, switching on the fields G in- 
creases all the moments of inertia, but particularly J7m, 
with respect to the case N . This causes a decrease in 
w^'ft , but the corresponding I^!^ does not change much, 
because Jm is larger. Switching on the fields T results in 
values of Jm between those obtained for the cases N and 
G. The critical frequency always becomes higher than for 
the G fields, and the resulting critical spin is always the 
highest among all the examined sets of time-odd fields. 
These variations in I^!^ are of the order of a few spin 
units. 

In ^^^La, where the HF chiral solutions were found 
in Section IIII Gl values of the HF critical frequency and 
spin, Lj^.fj and are slightly higher than the classical 
estimates, w^rft ^''^'^ -^c'u ■ This can be understood on the 
basis of results of the perturbative search for the HF chi- 
ral solutions, performed in Section Fill Gl As illustrated 
in Fig. 13 the perturbative ratio Im/^rm representing the 
moment of inertia with respect to the medium axis, J7m, 
slightly falls with rotational frequency, which, accord- 
ing to Eq. (|22|l . causes the rise of the critical frequency. 
When the time-odd fields are included, values of cj^?J. and 
^crH vary similarly to those of wjjj.f^ and /^'f^f . 

The HF method used in the present study does not take 
into account the pair correlations. In order to include 
pairing in self-consistent calculations one would have to 
apply the Hartree-Fock-Bogolyubov (HFB) method in 
the TAG regime with two quasiparticle states blocked. 
Present codes do not allow for such a study, and a sys- 
tematic investigation of pairing has to be left for future 
analyses. On the one hand, one can expect that pairing 
effects may facilitate the calculations by removing the 
sharp level crossings that occur in the HF method, but 
on the other hand, convergence in presence of blocked 
states may be more difficult to obtain. Supposing that 
the HFB TAG results can be inferred from the HFB PAG 
calculations through a determination of the parameters 
of the classical model, one could gain some information 
on pairing effects by analyzing the HFB PAG moments 
of inertia and alignments. 

The critical frequency represents the transition point 
between planar and chiral rotation. The notion of crit- 
ical frequency was first introduced in Ref. |22| |. where 
the expression (|22|) for its value was derived from the 



simple classical model presented in Section Fill El How- 
ever, the occurrence of a transition from planar to chi- 
ral regime in the structure of chiral bands was clear al- 
ready from earlier investigations. Both within the PRM 
and TAG, this effect was obtained numerically but left 
without comment. This transition is abrupt in the semi- 
classical cranking model, but may be rather smooth in 
real nuclei. This is because the angular momentum vec- 
tor oscillates about the planar equilibrium below ovyt, 
which corresponds to non- uniform classical rotations |34j , 
while above Wcritj it can still tunnel between the left and 
^ht chiral minima, which represents chiral vibrations 

Because of these reasons, mean-field methods can pro- 
vide quantitative description of the bands in question 
only in the low-spin, planar regime, where there is only 
one minimum. In the chiral region, the mean-field ap- 
proach does not take into account the interaction between 
the left and right minima, which are exactly degenerate 
in energy, and the experimentally observed energy split- 
ting between the chiral partners cannot be calculated. It 
is argued in the literature that the mean-field chiral solu- 
tion can be viewed as a kind of average of the two partner 
bands, and thus mean trends can be compared. One can 
also speculate about the value of the critical frequency or 
spin. Description of the transition region is an interest- 
ing topic for study invoking techniques beyond the mean 
field, like the Generating-Goordinate Method. 

Figure ^ gives a comparison of experimental and cal- 
culated energies in ^^^La. Full symbols denote the ex- 
perimental yrast (circles), Bl (squares), and B3 (dia- 
monds) bands, discussed in Section IIII Bl Open circles 
and black crosses represent the HF TAG planar and chi- 
ral solutions, respectively. Their classical counterparts 
are marked with dashed and solid lines. The HF results 
for the critical spin, I^J^ — 15.9 — 18.5?i, are rather high 
as compared to the spin range, in which the bands Bl and 
B3 are observed. Yet, the classical estimate, I^l^ = 9.2h, 
evaluated for the Total-Routhian-Surface (TRS) PAG re- 
sults is already below that range _22j. This means that 
the inclusion of pairing in the calculations may be im- 
portant for correct interpretation of the data. However, 
from the closeness of experimental spins to the possible 
values of Icru one can suppose that, whichever of the 
bands Bl and B3 could eventually be interpreted as the 
chiral partner of the yrast band, the concerned spin re- 
gion may actually represent the transition between pla- 
nar and chiral rotation. Although no HF chiral solutions 
were found in ^^°Gs, similar conclusions can be drawn for 
that isotone on the basis of the classical estimates of the 
critical spin. In case of i34pj. g^j^^ ^^^Pm it is not clear 
whether the oblate or triaxial solutions should be taken 
for comparison with the bands observed in those nuclei. 

At low spins in ^^^La, where the supposed chiral part- 
ners have not been observed, the yrast band is well repro- 
duced by the HF planar solutions, particularly with the 
time-odd fields included. This is consistent with the sup- 
posed planar character of rotation at low spins. Roughly 
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at the spin where the chiral partners commence to be 
visible, the yrast band significantly changes its behav- 
ior, which can be attributed to entering into the chiral 
regime. In this spin region the HF results agree semi- 
quantitatively with experimental energies. 



V. SUMMARY 

To conclude, the first Skyrme-Hartree-Fock (Skyrme- 
HF) calculations with the Tilted- Axis Cranking were per- 
formed in ^^°Cs, ^^^La, ^^''Pr, and -'^'^^Pm in the search 
for self-consistent solutions corresponding to nuclear chi- 
ral rotation. Only the configuration T^h\^^2 ^^11/2 ' earlier 
assigned to the observed candidate chiral bands in those 
isotones, was considered. Two Skyrme parametrizations, 
SLy4 and SkM*, were used. Terms depending on time- 
odd nucleonic densities were either kept or excluded from 
the Skyrme energy functional. 

From the Principal- Axis-Cranking analysis it was con- 
cluded that the system in question can be modeled by two 
gyroscopes, representing the valence particle and hole, 
with spins stiffly aligned with the short and long axes of 
a triaxial rigid rotor, which stands for the core. Such a 
model was analyzed in the classical framework. This led 
to an important conclusion that chiral rotation can exist 
only above a critical angular frequency, given by a simple 
expression. 

The HF solutions representing planar rotation were 
found in all the considered N = 75 isotones, and chiral 
solutions were obtained in ^^^La. These solutions pro- 
vide the first proof based on fully self-consistent methods 
that nuclear rotation can attain a chiral character. In all 
cases, the self-consistent solutions agree surprisingly well 
with the results of the classical model, which means that 
the model faithfully represents salient features of the ex- 
amined phenomenon. 

It was found that the time-odd densities in the energy 
functional have no qualitative influence on the results, 
and change mainly the moments of inertia. The HF val- 
ues of the critical frequency are rather high as compared 
to the spin range in which the candidate chiral bands 
were observed in ^'^^La. The HF energies agree satisfac- 
torily with experiment only in the low-spin parts of the 
bands, where the rotation is supposed to be planar. In 
the chiral regime, the mean field is unable to reproduce 
the data precisely, and the agreement is only qualitative. 
It seems, also, that the experimentally observed bands 
actually represent a transition from planar to chiral ro- 
tation. 

The criteria used so far when attributing the chirality- 
partnership interpretation to the experimentally found 
rotational bands that are based on the 'small energy 
splitting' argument are clearly unsatisfactory on a long 
run. Numerous superdeformed band studies (followed 
by the normal-deformation studies) have shown that dif- 
ferent intrinsic configurations and thus strictly speaking 
different-shape nuclei may manifest nearly identical ro- 



tational bands. It is therefore necessary, to provide the 
experimental evidence going beyond just the energy mea- 
surements, first of all the accompanying electromagnetic- 
transition information. This could help excluding the 
mistake of interpreting e.g. the shape coexistence phe- 
nomena in terms of chirality - without providing extra 
sufficient conditions. In this paper we were not able to 
propose any clear-cut necessary-and-sufficient condition 
criteria to attribute the chirality label to the experimen- 
tal bands either. However, we do believe that the con- 
cept of the critical frequency discussed in detail in this 
paper provides a useful tool in interpreting the experi- 
mental results. The very fact that in some nuclei the 
self-consistent HF calculations do provide chiral solutions 
is highly non-trivial and very encouraging message in this 
field of research. 
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APPENDIX A: ROTATIONAL PROPERTIES OF 
SINGLE-PARTICLE STATES IN A 
Dj-SYMMETRIC POTENTIAL 

In this Appendix, we discuss elementary rotational 
properties of s.p. eigenstates of a mean-field Hamilto- 
nian, which is symmetric with respect to the D2 group. 
The group D2 comprises the time-reversal operation, T, 
three signature operations, R^, Ry, Ry, which are rota- 
tions through 180° about the three Cartesian axes, and 
products of the time-reversal and signature operations, 
which are called T-signatures, cf. Ref. for more in- 
formation about this group in the context of mean-field 
calculations. In most cranking solutions corresponding 
to quadrupole deformation, the group is a symmetry 
group of the s.p. Hamiltonian, h, of Eq. For a sin- 
gle Kramers pair in a fixed potential, we investigate the 
response of the s.p. angular momenta to a cranking fre- 
quency applied in an arbitrary direction. Our conclusions 
are based only on symmetry arguments, and are thus in- 
dependent on the particular implementation of the mean 
field. 

Irrespective of spatial symmetries, whenever the s.p. 
Hamiltonian is invariant under the time reversal, its spec- 
trum exhibits the two-fold Kramers degeneracy. We con- 
sider a single Kramers pair, whose states are denoted as 
l/i) and l/i), where 

f |m> = s^Im) , TIm) = spl^) ; (Al) 
is an arbitrary phase factor and = —s^. 
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All information about the matrix elements of the 
angular-momentum operator, J, between the states 
and \p) can be represented in a convenient way in terms 
of the real alignment vector, , and the complex decou- 
pling vector, D^, of the state They are defined as 



{^l\J\^,) 



(m|J|m> 



(A2) 



Although the components of the decoupling vector 
change their phases when \fi) and |/2) change theirs, the 
relative phases of those components do not depend on 
the phase convention. Since the angular-momentum op- 
erator is odd under the time reversal, it can be easily 
verified that 



(A3) 



If the s.p. Hamiltonian is symmetric with respect to the 
D2 group, it is possible to chose the states forming the 
Kramers pair as eigenstates of any of the three signature 
operators, Ri, where i = x,y, z, but only one at a time, 
because the signature operators do not commute among 
themselves, i.e., for i ^ j 



RiRj — 



k—x.y.z 



Rk 



(A4) 



This results, respectively, in three formally different 
pairs, (I Mi), I Mi)), which are just three different bases 
in the same two-dimensional eigenspace of h. 

We choose states |/ii) so that they correspond to eigen- 
values —i under the action of Ri, while the eigenvalues 
of \p,i) are +i. Multiplication rules (|A4|) allow to easily 
express eigenstates \nx), iMa), I My), ^-^d I Mi/) through lin- 
ear combinations of eigenstates \fiz) and \fiz)- By fixing 
the relative phase between states j^z) and j/l^) we obtain 
the following expressions: 



tJ-x) = 




+ \fiz)) , 


(A5) 


fix) = 




- \fiz)) , 


(A6) 


\tJ-y) = 


/|(lA^.>H 


-■i\P-z)) , 


(A7) 


\fiv) = 




> + , 


(A8) 



where ^/i = exp(i7r/4) and \/— I = exp(— i7r/4). These 
formulae allow to write J^" , D^'' , J^^ , and D^'" in terms 
of J^^ and D"', i.e.. 



J''^- = Re D^'' , 

= J^^ - Ini 
J^" = -Im D^"" , 



iRe D^"' 



(A9) 
(AlO) 
(All) 
(A12) 



The fact that j^i) and j^i) are eigenstates of Ri, to- 
gether with the transformation rules of the components, 



Jj, of the angular momentum operator under the three 
signatures. 



Rf j j Ri 



+jj for j = i 
-Jj for j ^ i 



(A13) 



induces limitations on the components, Jj" and Dj% of 
the alignment and decoupling vectors. Namely, 



non-zero for j — i 

for j ^ i 

for j = i 

non-zero for j ^ i 



(A14) 
(A15) 



In other words, .7^' is confined to the axis i and £)'^' to 
the plane perpendicular to that axis. 

From these relations and for the eigenstates defined as 
in Eqs. ljA5|l - ljA8|l . it follows that all the quantities Jj"' 
and Dj^ can be expressed through the three "diagonal" 
components, J^^ , of the alignment vector, 

J^- = ( J,^- ,0,0), D^-^ (0, , -iJ^' ) , (A16) 
J^" (0, , 0) , £>^« = (-i ,0,Jz') , (A17) 
J^^ = (0,0, J^'^ ) , D"^ ^ ( J^== , -iJ^;y , 0) . (A18) 

Thus, all the information about the angular-momentum 
matrix elements within a Kramers pair in the spectrum 
of a D|^-symmetric s.p. Hamiltonian is contained in three 
real numbers, J^=" , Jy^ , and J^' . 

The symmetry group D2 itself does not impose any 
conditions on the "diagonal" components, J/*'. However, 
these values can be further restricted if some other sym- 
metry is present. For example, if h is axially symmetric, 
say with respect to the z axis, then the states |mz), Ifiz) 
are eigenstates of , which leads to quantization of J^^ . 
In fact, J^' = +1/2, —3/2, because Rz — exp{—inJz), 
while in the adopted convention Rz\lJ-z) — ~*IMz)- For 
states \fix) and \^iy), defined by HA5|I and HA7|I. one eas- 
ily finds 



J^^ = iRe(Mz|i+ + J-IMz 



J^y = -Re(/i,|J+- J_|/2, 



(A19) 
(A20) 



where J+ — + iJy and J_ 



Jx 



iJy are the ladder 



operators, that increment and decrement the magnetic 
quantum number, J^, of an eigenstate, \Jz), of Jz, 



J+\Jz) 



\Jz + 1) , 



J-\Jz 



\Jz - 1) 



(A21) 



One can see, therefore, that the matrix elements in (|A19|) 
and (|A20I) can be non-zero only if j/i^) and \p.z) differ 
in Jz by one, that is if J^' = 1/2. In such a case, 
(/iz|j-|Mz) — 0, and J^"" — Jy^ . These results can be 
summarized as 



(J^ ,,Jy ,,Jz ) — 



(J^^^, J^^, jr") for = 1/2 



(0,0,J[fi 



for J,| 



3/2,... 

(A22) 
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The parameter 



25t deserve 



(A23) 



is not restricted by the above kinematic conditions. It 
is related to the standard decoupling parameter, a — 
— 2J^^, considered by Bohr and Mottelson They 
take such phases for the states l/x^) and \jlz} that the 
T-signature-j/, Ry = TRy, is the complex conjugation in 
the basis formed by these states. Since J+ is even under 
Ry , in that convention the matrix element in Eq. (|A23|) 
is real. 

We now consider the TAG for a single Kramers pair 
in a Dl'-symmetric potential. We make two simplify- 
ing assumptions. First, that the Hamiltonian h in the 
Routhian h' of Eq. does not change with rotational 
frequency (non-selfconsistent cranking). Second, that the 
considered Kramers pair has no coupling to other eigen- 
states of h through the angular-momentum operator (iso- 
lated pair). The TAG under such conditions becomes 
a two-dimensional diagonalization problem, that can be 
solved analytically. 

We use the basis of states \fiz) and Iftz)- For a degener- 
ate Kramers pair, h reduces to its eigenvalue, e. Matrix 
elements of the angular-momentum operator are defined 
by Eq. (|A18p . Altogether, matrix of the s.p. Routhian 
© takes the form 



h' 



h — ujj 



(A24) 



" e 0' 




e 





ZUJy J y 



It is easy to verify that the two eigenstates of this 
Routhian have opposite-sign mean angular-momentum 
vectors J, whose components read 



J^=±- 



(C.2(J^)2+^2(JM„)2+^2(JM.)2) 



1/2 



(A25) 



Equation (|A25|I constitutes the central point of discus- 
sion, because it defines the sought response of the s.p. 
angular momenta to rotation under the assumed condi- 
tions. Note that the dependence of Ji on uJi is non- linear 
in the general case. 

Values of alignments HA25P are undefined if, and only 
if, all the products ' vanish - in particular when lu — 
0. In such a case, the Routhian (|A24p is proportional to 
unity, and the mean angular momenta of its eigenstates 
depend on their (arbitrary) unitary mixing. 



Two extreme cases of the dependence 
particular attention. 



• If J^^ = " = J|'- = g, then Ji ~ guji/w, and J 
always orients itself along oj, already for infinitesi- 
mal oj. This can be dubbed soft alignment. 



• When only one of the parameters J^=" , Jy^ , J^'' 
is non-zero, say then Ji — J^'Sij, and J is 
independent of u). We call this stiff alignment on 
the j-th axis. 



In axial nuclei, precisely one two-fold degenerate sub- 
state of each deformation-split j-shell has jj^^^ =1/2 and 

J^^ 7^ 0, which represents the soft alignment. Accord- 
ing to Eq. HA22|) . all other necessarily have a vanish- 
ing decoupling parameter, and are thus stiffly aligned 
with the symmetry axis. For prolate shapes, the lowest- 
energy substate has j!^^^ = 1/2, and is soft, while for 
oblate shapes it is the highest substate. In triaxial nu- 
clei, values of the parameters Jf', where i = s,m,l cor- 
responds to the short, medium, and long principal axes, 
are equal to the s.p. alignments obtained from the one- 
dimensional cranking about the three axes. Indeed, for 
cranking about the axis i, the s.p. states are eigenstates 
of Ri. For example, from the results of Section llll Dl one 
can see that for the lowest /111/2 substates of a triaxial 
nucleus only is non-zero, while for the highest sub- 
states only Jf''' does not vanish. These alignments are 
thus stiff. Note that there are no states with stiff align- 
ment on the medium axis. The response to rotation of 
the middle /111/2 substates is soft, because all the three 
parameters, , J^" , J^^ , are non-zero. 

In realistic cranking calculations the symmetry argu- 
ments discussed here interplay with the fact that there is 
angular-momentum coupling between different Kramers 
pairs and that the mean field does change with rotational 
frequency. In the results of the present paper, however, 
the change of deformation induced by rotation is negligi- 
ble. The angular-momentum coupling of the lowest and 
highest /111/2 substates to other s.p. states is rather weak, 
what can be seen from the small curvature of their one- 
dimensional Routhians in Fig. [3 The stiff character of 
their alignments is fully confirmed in our self-consistent 
calculations, as discussed in Section IlII Dl and further in 
the paper. Investigation on how suitable the notion of 
soft and stiff alignment is in other physical cases remains 
a subject for further research. 
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